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Abstract. We briefly review^ the overlap formalism for chiral gauge theo- 
ries, the overlap Dirac operator for massless fermions and its connection to 
domain wall fermions. We describe properties of the overlap Dirac operator, 
and methods to implement it numerically. Finally, we give some examples 
of quenched calculations of chiral symmetry breaking and topology with 
overlap fermions. 



1. Overlap formula for the chiral determinant 

In the overlap formalism [1], the chiral determinant is obtained by em- 
bedding the Weyl fermion inside a Dirac fermion through a many-body 
problem. Let be two many-body Hamiltonians 

= -(at 4)H^(«; 

^Talk given by Urs M. Heller at the workshop "Lattice fermions and structure of the 
vacuum", October 5-9, 1999,xi Dubna, Russia. 
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Let |0±) be the ground states of obtained by filhng all the positive 
energy states of H='=. Then 

detC(^) ^ (0- |0+)^ (2) 

"Proof": |0— ) is obtained by filling all the positive energy states of 75, 
and |0+) by filling all the positive energy states of H+. They are of the 
form: 

IO-): (J), a.d |0+): J-((^rT^'^^)„J. (3, 

where C^Cu^ = fJ'f.Uj. and is the normalization. 

The overlap formula is formal and needs to be regulated. It is valid 
only in the limit m — 00 and one should think of m as a pre-regulator. 
The formula is strictly valid only for ratios of determinants since there is a 
gauge field independent normalization in the formula. 

H"*" need not have an equal number of positive and negative energy 
states and this happens for topologically non-trivial gauge fields, where the 
difference between the number of negative and positive energy states of H"*" 
is 2Q. Then detC(^) = 0! Furthermore, 

(0-K\...ayO+) or (0-|ai,...a,|^||0+) (4) 

will be non-zero, for Q > or Q < 0, respectively, resulting in fermion 
number violation. Potenital anomalies reside in the phase of |0+). We will 
be concerned only with vector gauge theories, where only |(0— |0-|-)p enters 
and the anomaly is trivially cancelled in this case. 



2. Lattice regularization 

On the lattice H~ — = 75 remains unchanged, while 

H+ ^ H+ ^ H^(m) = j,DU-m) = ( ^%\~^'^ -B^^lm) 

where C{U) is the naive lattice discretization of 0(^4) and B(C/) is the 
standard Wilson term (with r = 1). 

|0— ) is still as in (3). Let ^ ~ 5) unitary matrix that 

diagonalizes Hj, with the first and second "block-column" spanning the 
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subspaces of positive and negative eigenvalues, respectively. Then, for a 
vector theory, 

1(0- |0+)p = detadeta^ (6) 
This can be obtained as the determinant of the overlap Dirac operator [2] 

1 



l+75e(Hj) 



where e{x) denotes the sign function. To see this consider, 
1 



a (5 
7 5 



+ 75 



a —13 
7 —5 



a 
5 



(7) 



(8) 



Since V is unitary we have det V = det 5/ det and hence we obtain 

det Dov = det a det . (9) 
The overlap Dirac operator can be generalized to the massive case 



^ov(/x) = ^ [l + /X + (1 - /x)75e(H+)] 

where — 1 < // < 1 is related to the fermion mass by [3] 

mf = Z-^lx{l + 0{a^)). 
The propagator for external fermions is given by 

D-\l^) = {l-li)-^ \d-Xii)-i\, 



(10) 



(11) 



(12) 



i.e. it has a contact term subtracted, which makes the massless propagator 
chiral: {Z)-i(0),75} = 0. 

A massless vector gauge theory can also be obtained from domain wall 
fermions [4], where an extra, fifth dimension, of infinite extent is introduced. 
In the version of ref. [5], one can show [6] that the physical (light) fermions 
contribute log det .Ddw to the effective action with the 4-d action 



-Ddw 



1 + + (1 - //)75 tanh 



logT 



(13) 



where T is the transfer matrix in the extra dimension and Lg its size. As 
long as log T 7^ we obtain in the limit as — > oo 



i^DW ^ - [1 + /X + (1 - /x)75e(- log T)] 



(14) 
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This is just the massive overlap Dirac operator up to the replacement — > 
— log T. It is easy to see that in the limit ag — 0, where ag is the lattice 
spacing in the extra dimension (set to 1 above), one obtains — logT = 
H^(l + 0(a,)). 

3. Some properties of the overlap Dirac operator 

In many cases it is more convenient to use the hermitian version of the 
overlap Dirac operator (10): 

Hoiii) = 75^ov(/x) = ^ [(1 + /^)75 + (1 - tJ^)e{H^)] . (15) 

The massless version satisfies, 

{i/„(0),75} = 2//2(0). (16) 

It follows that [ii'o(0),75] = 0, i.e. the eigenvectors of i^o(O) can be chosen 
as chiral. Since 

Hl{i,) = {I - ^?)Hl{Q) + ix-" (17) 

this holds also for the massive case. 

The only eigenvalues of Ho{G) with chiral eigenvectors are and ±1. 
Each eigenvalue < < 1 of H^{0) is then doubly degenerate with 
opposite chirality eigenvectors. In this basis Ho{iJ,) and Dov(m) block 
diagonal with 2x2 blocks, e.g 

( (i-^)A2 + ;. (i-^)Ayr^\ 



where 



75= (J (19) 



For a gauge field with topological charge Q 7^ 0, there arc, in addition, 
\Q\ exact zero modes with chirality sign((5), paired with eigenvectors of 
opposite chirality and eigenvalue 1. These are also eigenvectors of Ho{fi) 
and Dov(m)- 

-Cov ( At) zero sector : (oi) °^ (o/i 

(20) 

depending on the sign of Q. 

We remark that from eigenvalues/ vectors of Hg{0) those of both Ho{fi) 
and DoY{fJ,) are easily obtained. There is no need for a non-hermitian eigen- 
value/vector solver! For example, the Ritz algorithm [7] will do just fine. 
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4. Small eigenvalues and the chiral condensate 

In the chiral eigenbasis of Hg{0) the external propagator takes the block 
diagonal form with 2x2 blocks 















2) 


or 


(o 


°) 



^ ■ A2(l-^2) + ^2 Ax/r^ m(1 - j ' ^^^^ 

and, in topologically non-trivial background fields the \Q\ additional blocks, 
depending on the sign of Q, 

\ / n n \ 

(22) 

We thus find in a fixed gauge field background 

and averaged over gauge fields wc get the condensate. It is dominated by 
the small (non-zero) eigenvalues and in the thermodynamic limit, where 
the first term vanishes, it is given by the density of eigenvalues at zero, 
p(0+). 

With our normalizations we find for all chiral vectors \b) 

ti{b\[r5D-Hl^)f\b) = {b\D-\ij)\b) V6 with 75|6) = ±|6). (24) 

This ensures the relation fiXn = 2(V'V') ^or every configuration, and, in fact, 
for every chiral random source used in a stochastic estimation of condensate 
and chiral susceptibility x-k- For such stochastic estimates, we always work 
in the chiral sector with no zero-modes. 



5. Implementations of the overlap Dirac operator 

In practice, we only need the application of D(ii) on a vector, D{fj,)ip, 
and therefore only the sign function applied to a vector, e{Hw)il^- Since we 
need the sign function of an operator (a large sparse matrix) this is still a 
formidable task. 

Methods proposed for this computation are: 

— A Chebyshev approximation of e(x) = over some interval [5, 1] [8]. 
For small 5 a large number of terms are needed. 

— A fractional inverse method using Gegenbauer polynomials for [9] . 
This has a poor convergence since these polynomials are not optimal 
in the Krylov space. 
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Figure 1. Plots of the optimal rational function approximation to e{x) for various order 
polynomials. 



— Use a Lanczos based method to compute based on the sequence 

generated for the computation of ^ [10]. 

— Use a rational polynomial approximation for e(x) which can then be 
rewritten as a sum over poles: 

The application of x ^ ^iHw)ip can be done by the simultaneous 
solution of the shifted linear systems [11] 

{Hl + bk)4>k = ''P, X = Hu,{cQij + Y,Ck(t>k)- (26) 

k 

One such approximation, based on the polar decomposition [12], was 
introduced in this context by Neuberger [13]. We use optimal ratio- 
nal polynomials [14]. The accuracy of this approximation is shown in 
Fig. 1. 

We note that in all methods listed above, one can enforce the accuracy 
of the approximation of e{x) for small x by projecting out the lowest few 
eigenvectors of Hy^ and adding their correct contribution exactly. 



i=l 4=1 



(27) 

To invert D'^ D for overlap fermions, we have, generically, an outer 
CG method (a 4-d Krylov space search) and an independent inner search 
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Figure 2. The approach of p(0; m) to the continuum hmit in the quenched theory. 



method for €{Hyj)tp - maybe CG again. For domain wall fermions, on the 
other hand, a 5-d Krylov space search method is used. It may pay off to 
try to combine inner and outer CGs for overlap fermions by reformulating 
them into a 5-d problem [15, 16]. 

6. Main problem for Overlap and Domain Wall fermions 

For topology to change, we must create dislocations. These produce small 
modes which force the spectral gap of if^(m) to be closed. The density 
of zero eigenvalues of Hyj{m), p{0;m), is non-zero in the quenched case, 
but rapidly decreasing with decreasing coupling [17]. Very roughly, we find 
p{0;m)/a^^^ ~ e~'^'^ as shown in Fig. 2. 

The existence of small eigenvalues hampers the approximation accuracy 
and convergence properties of implementations of e{Hyj). Eigenvector pro- 
jection both increases the accuracy of the approximation and decreases the 
condition number, e.g. of the inner CG. 

The existence of small eigenvalues has implications also for domain wall 
fermions. One can show that the spectrum of — log T{m) of Eq. (13) around 
zero is the same as the spectrum of Hyj{m) [1]. While the small eigenval- 
ues of — log r(m) don't appear to cause algorithmic problems for domain 
wall fermions, they can induce rather strong Lg dependence of physical 
quantities, and causing hence the need for large Lg. 



o m = 1.2 
o m = 1.5 
X m = 1.7 
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Figure 3. Examples for the microscopic spectral density for all three ensembles. For 
overlap fermions one can probe different topological sectors (rightmost panel). 



7. The Dirac spectrum, chiral condensate and chiral Random 
Matrix Theory 

Up to a scale, given by the infinite volume chiral condensate E = {il^ip)-, 
RMT predicts that the rescaled density of eigenvalues 

is universal, dependent only on the symmetry properties, number of dynam- 
ical flavors, and the number of exact zero modes (the topological sector), 
but not the form of the potential in the random matrix theory, or low 
energy effective Lagrangian [18]. There are three classes of random ma- 
trices determined by their symmetry properties: orthogonal, unitary, and 
symplectic. 

In Fig. 3 we show examples of the microscopic spectral density for all 
three ensembles and compare to the analytic predictions from RMT. With 
overlap fermions we can probe topologically non-trivial sectors. 

Similarly, there are predictions in each ensemble and topological sector 
for the distribution of the lowest eigenvalue. Examples for the quenched 
theory with overlap fermions are shown in Fig. 4. The S's extracted from 
fits in different v sectors for each ensemble are consistent [19]. 

RMT also gives predictions for the finite mass and volume dependence 
of the chiral condensate in the small mass large volume regime. 



Hdz-^ +^ (29) 



S 

with u = /xSy. Particularly interesting is the behavior at small u: 
sGUE(^)/S^-^logu, sGOE(^)/s^l(;r-u), Sgf (u)/S ~ u (30) 
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It is very sensitive to the lowest eigenvalues. In quenched QCD, surprisingly, 
Sg'^^(M)/S does not vanish at u = in this microscopic limit. Our data, 
shown in Fig. 5 follow the predictions well [20]. Once again, with overlap 
fermions we can probe topologically non-trivial sectors. The finite volume 
corrections are quite large for overlap fermions: E is about a factor 7 smaller 
than in the staggered case. This implies that for overlap fermions larger 
volumes are needed to see the microscopic regime. 

8. Small eigenvalue distribution in quenched QCD above Tc 

We have studied the small eigenvalue distribution of the Dirac operator 
in the deconfined phase of quenched QCD. Sample distributions of small 

(non-zero) eigenvalues are shown in Fig. 6 [21]. 

For overlap fermions, we see the lower end of the bulk of the distribution, 
then a dip, or even a gap, and then again small eigenvalues, below about 
0.05. We focus on the small modes, A < 0.05. Our findings are summarized 
in Tables 1 and 2. We see that both {n)/V and {Q'^)/V seem to remain 
finite and non-zero in the large volume limit for fixed /3, but they drop 
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Figure 6. Low lying eigenvalue distributions in quenched QCD at finite temperature for 
staggered and overlap fermions. 



quickly as /3, and hence the temperature, is increased. 
Looking in more detailed at the small modes we find 

— Their numbern is roughly Poisson distributed, P(n, (n)) = (n)"e~^"^/ra!. 
Average and variance are approximately equal. 
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TABLE 1. SU(3) data: n = n^+ri- with n± the number 
of zero and small non-zero eigenvalues with chirality ±. 
Q = n+—n- is the topological charge. cr„ is the variance 
of n. The volume normalizations for n and are per 
spatial 8^ volume. 



volume 


8^ X 4 
5.75 5.85 


12^ X 4 
5.75 


16^ X 4 
5.71 5.75 5.85 


{n)/V 

{Q')/v 

(n) 1 cr„ 


0.32 0.06 
0.31 0.07 
1.09 0.90 


0.28 
0.28 
0.92 


0.63 0.30 0.05 

0.64 0.33 0.05 
1.15 1.03 0.83 



TABLE 2. SU(2) data 



volume 

/3 



2.3 



X 4 
2.4 



le-" 

2.4 



X 4 
2.5 



(n)/l/ 
{n)/a„ 




0.25 0.05 
0.25 0.05 
0.93 0.99 



— For fixed n, n+ and n_ are roughly binomially distributed. 

These observations are consistent with interpreting the small modes to be 
due to a dilute gas of instantons and anti-instantons, with n+ and n_ their 
numbers, n — \Q\ of the would-be zero modes mix due to their overlapping 
and get small eigenvalues, while \Q\ exact zero modes remain. 

At finite temperature, instantons fall off exponentially, and so do the 
fermionic zero modes associated with them. We consider a toy model of ran- 
domly (Poisson and binomially) distributed instantons and anti-instantons, 
inducing interactions of the form hoe~'''^'''^^^^ between the would-be zero 
modes of every instanton - anti-instanton pair (i, j) with separation d{i,j). 
Like sign pairs are assumed to have no interactions. This toy model repro- 
duces all qualitative features of the small eigenvalue distributions well for 
D 2, corresponding to D Ki 1/{2T) [21]. 

9. Conclusions 

The overlap Dirac operator has the same chiral symmetries as continuum 
fermions. It has exact zero modes in topologically non-trivial gauge fields. 
It is therefore well suited for a study of the interplay of topology, with its 
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associated exact zero modes, and chiral symmetry breaking, determined by 
the density of small eigenvalues. 

In the range of its validity the predictions of chiral random matrix theory 
are well followed and confirmed by overlap fermions, including the depen- 
dence on topology, given by the number of exact zero modes. 

A study of the small eigenvalues in quenched QCD above the deconfining 
transition temperature, Tc, shows that topology, manifested by exact zero 
modes, persists. Furthermore, a finite density of small eigenvalues persists, 
and their properties are well described by attributing them to the would-be 
zero modes of a random dilute gas of instantons and anti-instantons. 
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